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ABSTRACT 

All but three (M87, BL Lac and 3C 279) extragalactic sources detected so far at very high 
energy (VHE) 7-rays belong to the class of high-frequency peaked BL Lac (HBL) objects. This 
suggested to us a systematic scan of candidate sources with the MAGIC telescope, based on 
the compilation of X-ray blazars by Donato et al. (2001). The observations took place from 
December 2004 to March 2006 and cover sources on the northern sky visible under small zenith 
distances zd < 30° at culmination, constraining the declination to values between —2° and 
+58°. The sensitivity of the search was planned for detecting X-ray bright (^(lkeV) > 2/iJy) 
sources emitting at least the same energy flux at 200 GeV as at 1 keV. In order to avoid strong 
7-ray attenuation close to the energy threshold, the redshift of the sources was constrained to 
values z < 0.3. Of the fourteen sources observed, 1ES 1218+304 (for the first time at very high 
energies) and 1ES 2344+514 (strong detection in a low flux state) have been detected in addition 
to the known bright TeV blazars Mrk 421 and Mrk 501. A marginal excess of 3.5ct from the 
position of 1ES 1011+496 was observed and has been confirmed as a source of VHE 7-rays by 
a second MAGIC observation triggered by a high optical state ( Albert et al. 2007ef ). For the 
remaining sources, we present here the 99% confidence level upper limits (u.l.) on the integral 
flux above ~ 200 GeV. We characterize the sample of HBLs (including all HBLs detected at VHE 
so far) by looking for correlations between their multi-frequency spectral indices determined 
from simultaneous optical, archival X-ray, and radio luminosities, finding that the VHE emitting 
HBLs do not seem to constitute a unique subclass. The absorption corrected 7-ray luminosities 
at 200 GeV of the HBLs are generally not higher than their X-ray luminosities at 1 keV. 



1. Introduction 

Blazars belong to the most extreme objects 
in astronomy. Dominated by a non-thermal con- 
tinuum spectrum, covering up to twenty decades 
in energy, they show variab ility on time scales 



length from the nucleus indicate in situ parti- 



cle acceleration (jBiretta et al.j|1991l) . occuring at 



of years down to min utes ([Albert et al. 2007f; 



Aharonian et al.l 12007a ) , and apparent luminosi- 
ties exceeding 10 49 ergs _1 . Morphologically, 
blazars show strongly collimated jets extending 
from scales not much larger than t he event hori- 



zon o f a supermassive black hole (jBiretta et al 



12002! ) up to megaparsec scales. Superluminal mo- 
tion of knots in the radio jets indica tes relativistic 
bulk motion (Ghiscllini et al. I ll993l) . X-ray knots 
at distances of more than the radiative cooling 
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traveling and stationary shocks i n the jet. Accord- 
ing t o the unified scheme (e.g. Illrrv fe Padovani 
1995 ), blazars are accreting supermassive black 



holes expelling a relativistic plasma jet at a small 
angle between the jet axis and the line of sight, 
with strong boosting of the observed emission due 
to relativistic bulk motion of the emitting plasma. 
BL Lacertae objects differ from the generally more 
luminous quasars by showing only faint or even ab- 
sent emission lines, the absence of thermal big blue 
bump emission, and by not showing the otherwise 
typical luminosity evolution. 

The Spectral Energy Distribution (SED) of 
blazars shows two pronounced peaks, the first be- 
tween IR and hard X-rays, which is commonly 
believed to be synchrotron radiation of highly 
relativistic electrons, and the second one at 7- 
rays. Depending on the location of the first 
peak, BL Lac objects are further divided into low- 
frequency peaked BL Lacs (LBL, IR to optical) 
and high- frequency peaked BL Lacs (HBL, UV to 
X-rays) (jGiommi fc Padovanil ll994). The second 
peak at high energies can be explained by in- 
verse Compton scattering of low energy photons, 
produced as synchrotron radiation by the same 
populatio n of electrons (Synch rotron Self Comp- 
ton, SSC. lMaraschi et al.l (|l992r i). or from ambient 
thermal photon fields, whi ch could enter directly 
into the emission region (jDermer fc Schlickeisei 



19931 ) or by scattering on material surrounding 



the jet (jSikora et al.l 11994 ). The origin could 



also be due to hadronic processes associated with 
proton and ion acceleration which leads to elec- 
tromagneti c cascades and p roton synchrotron 
radiation dMannheim Il993t Aharonian 2000; 
Muecke fc Protherod 2001 ). We cautiously re- 



mark, that the SED is probed at a sufficient level 
of sensitivity only in a limited range, there are 
still large gaps, in particular between 50 keV and 
100 MeV where further peaks could show up. 

In December 2004, when the regular observa- 
tions with the Major Atmospheric Gamma-ray 
Imaging Cherenkov (MAGIC) telescope started, 
the number of known VHE blazars was six, all of 
them X-ray bright HBL objects. At the time of 
writing, the number has increased to 193, includ- 



1 For an up-to-date list of VHE blazars, see 
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ing on e LBL object (BL Lacertae, lAlbert et al 
2007 dl), one Flat Spe ctrum Radio Quasar (3C 279 



Teshima et al.l 20071) and the giant radio galaxy 



M87 (jAharonian et alj|2006bl ). 

The detection of VHE 7-rays from cosmologi- 
cal distances is made difficult, due to absorption 
of 7-rays by photon-photon interactions with low 
energy photons from the evolving metagalactic ra- 
diation field (MRF). In the 100 GeV to 10 TeV 
range, far-infrared to optical photons are most im- 
portant for the attenuation. It has been realized 
that this leads to a relation between the 7-ray 
cutoff energy and the sou rce redshift known as 
the Fazio-Stecker relation ( Fazio fc Ste ckcr 1979; 
iKneiske et al.|[20o3 ). The fact that the so-far de- 
tected VHE sources have much lower redshifts 
compared with the EGRET 7-ray sources is in line 
with the expected effect of 7-ray attenuation, al- 
though the lack of curvature in t he observed spec- 
tra is a source of serious doubts ( Aharonian et al 



2006aD . 

Due to the small field-of-view of an Imaging 
Air Cherenkov Telescope (IACT) and the limited 
duty cycle of ~ 1000 firs per year, promising can- 
ditates for VHE emission have to be selected care- 
fully. All established TeV sources are bright X-ray 
sources, most of them with comparable luminosi- 
ties in both regimes, which renders a systematic 
scan of the X-ray brightest HBL objects a reason- 
able approach. 

Here we report on the results of such an ap- 
proach pursued with the MAGIC telescope for a 
sample of 14 HBLs. In Sect. 2 the selection criteria 
for the sample are discussed, while the description 
of the observations can be found in Sect. 3. The 
data analysis technique is described in Sect. 4, and 
the analysis results are summarized in Sect. 5. A 
brief explanation of how the SEDs were obtained 
from the data (using archival radio and X-ray data 
as well as simultaneous optical data) can be found 
in Sect. 6. Finally, the resulting properties of the 
SED of HBLs and inferences on their luminosity 
function are discussed in Sect. 7. 

2. HBL sample 



We used the compilation from iDonato et al 



(|200l[ ). which provides 421 X-ray fluxes with spec- 



tral information for 268 blazars (136 of them HBL 
objects) together with average radio (at 5 GHz) 
and optical (V-band) fluxes. The selection cri- 
teria were (i) redshift z < 0.3, (ii) X-ray flux 
F x (l keV) > 2 juJy, and (iii) zenith distance (zd < 
30° during culmination, a total of fifteen sources 
remained after cuts. 

The selection was made to avoid strong 7-ray 
attenuation at the energy threshold. At z = 0.3, 
the expected cut-off energy is still above 200 GeV, 
where MAGIC has its highest sensitivity. As the 
energy threshold increases with the zenith dis- 
tance, all observations were carried out below 
40°, where the analysis energy threshold is around 
200 GeV. As most of the established TeV sources 
show comparable luminosities in X-rays and in 7- 
rays, only the X-ray brightest HBLs were selected, 
leading to a cut at 2 /iJy. Assuming the same lu- 
minosity at ~ 200 GeV, it corresponds to ~ 7% of 
the flux of the Crab Nebula, which would be de- 
tectable for MAGIC within 15 hrs. 

The goal was to observe them for at least 15 hrs 
in order to establish new VHE sources and to put 
constraints on the SED of HBLs in a systematic 
fashion. The complete set is listed in TablcQ] 

3. Observations 

The MAGIC telescope is a single dish IACT, lo- 
cated on the Canary island of La Palma (28.8° N, 
17.8° W, 2200 m a.s.L). A 17 m diameter tes- 
sellated parabolic mirror with a total surface of 
234 m 2 , mounted on a light-weight space frame 
made from carbon fiber reinforced plastic tubes, 
focuses Cherenkov light from air showers, initi- 
ated by 7-rays or charged cosmic rays, onto a 
576-pixel photomultiplier camera with a field-of- 
view of 3.5°. The analogue signals are trans- 
ported via optical fibers to the trigger electronics 
and each channel is read out by a 300 Msample/s 
FADC. Further details on th e tel escope can be 
found in Baixeras et al. I (|2004h and lCortina et all 
(|2005l ) . Note that the readout system has been up- 
graded to a 2GSamp le/s FADC in February 2007 



(jGoebel et al.ll2007h . A second telescope of the 



http: / /www.mppmu.mpg.de~rwagncr/sources/ 



same size for observations in stereo mode is cur- 
rently under construction. 

The observations took place from December 
2004 to March 2006 in moonless nights. The data 
are taken in different observation modes. If the 
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Table 1 
List of targets. 



Source 




RA 




dec 




z 


flux a 


r b 




season 


Time [h] c 


1ES 0120+340 


01 


23 


08 


9 


+34 20 50 





272 


4.34 


1 


93 


2005 


Aug-Sep 


14.9 


RX J0319.8+1845 


03 


19 


51 


8 


+ 18 45 35 





190 


1.76 d 


2 


07 


2004 


Dec-2005 Feb 


6.9 
























2005 


Sep-2006 Jan 


4.7 


1ES 0323+022 


03 


26 


14 





+02 25 15 





147 


3.24 


2 


46 


2005 


Sep-Dec 


11.4 


1ES 0414+009 


04 


16 


53 




+01 04 54 





287 


5.00 


2 


49 


2005 


Dec-2006 Jan 


17.8 


1ES 0806+524 


08 


09 


49 


2 


+52 18 58 





138 


4.91 


2 


93 


2005 


Oct-Dec 


17.5 


1ES 0927+500 


09 


30 


37 


6 


+49 50 24 





188 


4.00 


1 


88 


2005 


Dec-2006 Mar 


16.1 


1ES 1011+496 


10 


15 


04 


2 


+49 26 01 





212 c 


2.15 


2 


49 


2006 


Mar-Apr 


14.5 


Mrk 421 


11 


04 


27 


3 


+38 12 31.8 





030 


39.4 


2 


96 


2004 


Nov-2005 Mar 


25.6 f 


1ES 1218+304 


12 


21 


21 


9 


+30 10 37 





182 


8.78 


2 


34 


2005 


Jan 


8.2 s 
























2006 


Jan-Mar 


14.6 


RX J1417.9+2543 


14 


17 


56 


6 


+25 43 25 





237 


3.58 


2 


25 


2005 


Apr-Jun 


13.0 


1ES 1426+428 


14 


28 


32 


5 


+42 40 25 





129 


7.63 


2 


09 


2005 


Mar-Dec 


6.1 


Mrk 501 


16 


53 


52 


2 


+39 45 36.6 





034 


20.9 


2 


25 


2005 


May- Jul 


32.2 h 


RX J1725.0+1152 


17 


25 


04 


4 


+ 11 52 16 


> 


0.17' 


3.60 


2 


65 


2005 


Apr 


5.3 


1ES 1727+502 


17 


28 


18 


6 


+50 13 11 





055 


3.68 


2 


61 


J 




j 


1ES 2344+514 


23 


47 


04 


9 


+51 42 18 





044 


4.98 


2 


18 


2005 


Aug-2006 Jan 


23. l k 



a F(l keV) [/tJy], average value from different measurements, taken from lDonato et alj (|200ll ). 
b r is the spectral index for the differential spectrum (dN/dE) at IkeV, assuming a power law 
c Effcctive observation time after quality selection 
d Two measurements are above 3 /xJy, one below 1 /iJy 

°The earlier reported redshift of 0.200 was recently revised by I Albert et al. ( 2007el) 
f Results published in lAlbert et all (|2007d) 



s Results published in I Albert et all (|2006bl ) 



h Results published in I Albert et all |2007J) 



'The earlier reported redshift of 0.018 was recently revised by a lower limit (jSbarufatti et al.l 200£ 
1 Proposed but not observed due to bad weather 
k Results published in I Albert etHI (|2007bh 
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telescope is pointing to the source (on- mode), the 
background has to be determined by so-called off- 
data, where the telescope points to a nearby sky 
region where no 7-ray source is expected. The 
off-data cover the same zd range with a similar 
night-sky background light intensity. The larger 
fraction of the source sample was observed in the 
so-called wobble mode, where the pointing of the 
telescope wobbles every 20 minutes between two 
symmetric sky locations with an angular distance 
of 0.4° to the source. The background in the signal 
region can be estimated from sky locations placed 
at the same distance from the camera center as 
the candidate source. 

Except 1ES 0927+500 and 1ES 0414+009, all 
objects were monitored by the KVA telescope 
(http://users.utu.fi/kani/lm/index.html) on La 
Palma in the optical R-band. None of the sources 
showed flaring activity in the optical during the 
MAGI C observations. Th e host galaxy corrected 
fluxes (jNilsson et al.ll2007l ). taken simultaneously, 
averaged over the time of the MAGIC observa- 
tions, are listed in TableO 

4. Data analysis technique 

The data are processed using the MAGIC Anal- 
ysis an d Reconstruction Software (MARS) ( Bret j 
2005ah . A descriptio n of the different analysis 
steps can be found in iGaug et all (T2 005) (includ- 
ing the calibration) and Bretd ( 2005bh . 

As the trigger rate strongly depends on the 
weather condition, only data with a rate above 
160 Hz are used to ensure a high data quality. 

The moments up to third order of the light dis- 
tribution are used to cha racterise eac h event by a 
set of image parameters (|Hillaslll985l) . For back- 
ground suppression, a SIZE-dependent pa rabolic 
cut i n WIDTH x LENGTH is applied (jRiegell 
2005). To reconstruct the origin of the shower 



tion of the shower images at the camera border 
is taken into account. These coefficients are de- 
termined separately for on-off and wobble data 
using simulated 7-showers, wh ich are produced 



in the camera plane, the DISP method is em- 
ployed (jLessard et al.l 120011 ) to estimate the dis- 
tance between the centre of gravity of the shower 
and its origin. The third moment determines 
the direction of the shower development. The 
constant coefficient £ from the parametrisation 
of DISP in the original approach is replaced in 
this analysis by £ + £1 • (LEAKAGE) 6 , LEAK- 
AGE being the fraction of light contained in the 
outermost camera pixels. Thereby the trunca- 



by CORSIKA, versio n 6.023 dHeck et alJ fl998: 



Maiumdar et al. 2005) for zd below 40° and en- 



ergies between lOGeV and 30TeV, following a 
power law with a spectral index -2.6. 

The cut coefficients for the background suppres- 
sion are optimised using Crab Nebula data, taken 
at similar zd. One set of cut coefficients is derived 
for data taken in on-off mode and one for wob- 
ble mode. The significance of a possible signal is 
determined from the distribution of the squared 
angular distance (9 2 ) between the shower origin 
and the source position. The signal region is de- 
termined as 9 < 0.23°, corresponding to slightly 
more than two times the 7 point-spread function 
of the MAGIC telescope. 

For observation in on-off mode, the off-data 
have to be scaled to match the on-data. This is 
done in the region 0.37° < 9 < 0.80°, where no 
bias from the source is expected. For wobble ob- 
servations three regions, located symmetrically on 
a ring around the camera centre with the same dis- 
tance from the centre as the source position, are 
defined as background regions. The scale factor is 
fixed to 1/3. 

For every source the sta tistical significa nce ac- 
cording to equation 17 from Li fc Ma ( 19831 ) is cal- 
culated. The upper limit on the excess rate is de- 
rived with a con fidence level ( c.l.) of 99%, using 
the method from lRolke et alJ (2005), which takes 
also the scaling factor of the background into ac- 
count. The upper limit for the excess rate is then 
compared to the excess rate of the Crab Nebula, 
which leads to an upper limit of the flux in units 
of the Crab Nebula flux above a certain energy 
threshold, assuming a Crab like spectrum. The en- 
ergy threshold is here defined as the energy where 
the differential distribution (dN/dE vs. E) of simu- 
lated 7-ray events, surviving all cuts, peaks. Note 
that this threshold depends also on the spectral 
shape. 

A large sample of data from the Crab Nebula 
in on-off as well as in wobble mode is used, spread 
over the entire observation campaign (see Tablc[2]). 
This analysis shows, that the excess rate of the 
Crab Nebula is correlated to the background rate 
(after 7-hadron separation). Therefore, depend- 
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Table 2 

Observations of the Crab Nebula used for the upper limit calculation 







Exp. 


zd 


~Ethres 


Excess 


Background 




Sign. 


Season 


Mode 


[hrs] 


degr. 


GeV 


[mill -1 ] 


[min -1 ] 


scale 


[a/Vh] 


2005 Oct - 2006 Mar 


w 


8.6 


23.0 


230 


6.43 


6.23 


0.33 


14.6 


2004 Dec - 2006 Mar 


on 


38.7 


16.2 


190 


8.80 


9.19 


0.92 


13.5 




Fig. 1. — Excess rate (excess) vs. background rate 
(bgd) for the Crab Nebula (on-off mode). A linear 
fit yields (0.95 ± 0.10) ■ bgd + 0.2 = excess. 



Fig. 2. — Excess rate (excess) vs. background 
rate (bgd) for the Crab Nebula (wobble mode). A 
linear fit yields (1.26 ± 0.3) • bgd — 1.5 = excess. 



ing on the background rate of the AGN, a refer- 
ence value for the excess rate of Crab has to be 
calculated. This can be understood when taking 
into account that even after quality selection the 
rates fluctuate up to 20%, depending on weather 
conditions. In Fig. Q] and Fig. [5] the rate of excess 
events vs. the background rate of different subsam- 
ples of the Crab Nebula is shown. A linear regres- 
sion to the on-off samples results in an acceptable 
fit (% 2 /d.o.f. = 26/14), showing a clear correla- 
tion. The fit for the wobble data is quite poor 
(% 2 /d.o.f. = 13.5/3). As a constant fit gives even 
a worse result and the fit on the on-data shows a 
correlation between background- and excess rate, 
the linear fit for the wobble data is also used to 
calculate the reference values for the comparison 
of the excess rates. The Crab units are converted 
into a flux of photons cm " 2 s -1 using the spect rum 
of the Crab Nebula from I Albert et all (|2007gh . 
The systemat ic error for the flux is estimated 
Albert et all (|2007g|) and discus- 



sion therein). For the u.l. determination there is 
also the uncertainty of the correct energy thresh- 
old (which depends on the source spectrum). 

5. Results of the MAGIC observations 

Within this observation program, VHE 7-rays 



were d iscovered from 1ES 1218+304 (jAlbert et al 
2006bh and 1ES 2344+515 was observed in a 



low flu x state with high significance (jAlbert et al 
2007bl ). Mrk 421 was observed for more than 
25 hrs in 2005. The results are discussed in detail 



in lAlbert et all (|2007d) . Mrk 501 was observed 
from May to July 2005 with more than 30 hrs, re- 
vealing a high precision lightcurve on a day-by-day 
basi s as well a two exce ptionally short-time flares 
(see lAlbert et al.l (|2007t ) for more details). For ten 



to be -30% (see 



sources of the sample, no significant signal is seen. 
The 2006 observation of 1ES 1218+304 results in 
a weak signal of 4.6 a (see Sect. 5.2). A slightly re- 
fined analysis of 1ES 1011+496 yields a hint of a 
signal with a 3.5 a significance (see Sect. 5.3). The 
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results are listed in Table[3] Observations of 1ES 
1727+502 are still pending. 

5.1. Upper limits 

The u.l.s are between 2.3% and 8.6% of the 
Crab Nebula flux. For a Crab-like spectrum the 
energy thresholds vary between (190±15) GeV and 
(230 ± 15) GeV, depending on the zd of the obser- 
vation. For the threshold calculation the exact 
zd distribution of every observation is taken into 
account. As the Crab spectrum at ^200 GeV is 
quite hard (spectral slope -2.26 for the differen- 
tial energy spectrum), the u.l.s are also calculated 
for an -3.0 power law spectrum, which represents 
quite well the average slope of all HBLs detected 
at VHE so far. In Table[5]the energy threshold as 
well as the flux u.l. at 200 GeV are given under the 
assumption of a -3.0 power law. 

5.2. 1ES 1218+304 

The source was observed in 2006 from January 
the 29th to March the 5th during 15 nights with in 
total 14.6 hrs. Figure[3] shows the distribution of 
the squared angular distance between the recon- 
structed shower origin of each event and the as- 
sumed source position. The vertical line indicates 
the signal region. The background is determined 
by three off-regions in the camera. The excess has 
a statistical significance of 4.6 a (see also Table[3|). 



trum with spectral slope of -3.0 as measured in 



2005 (| Albert et all l2006bl) . the energy threshold 




l,uu 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 
e 2 [degr 2 ] 

Fig. 3. — Distribution of the squared angular dis- 
tance with respect to the position of 1ES 1218+304 
(dots) and to three off- regions (scaled by 1/3, gray 
shaded area). 

Under the assumption of a power law spec- 



decreases with respect to a Crab like spectrum to 
190 GeV. The average integral flux above 180 GeV 
for the complete sample is F(> 180 GeV) = 
(1.48 + 0.48) x lO" 11 photons cm^s" 1 . 

The search for flux variability on time scales 
from days (~ 1 h observation time per night) to 
weeks (sample with several consecutive nights) 
yielded no significant variability. 

The integral flux above 180 GeV indicates a 
~ 30% lower flux than the one measured in 2005, 
even though they are consistent within their er- 
rors. The average optical flux in 2006 was ~ 20% 
lower than in 2005, which is already significant 
compared to the statistical error of ~ 2%. During 
the observations in 2006 the optical flux decreased 
continuously from (1.144 ± 0.036) mJy on the 3rd 
of February to (0.947 ± 0.038) mJy on the 7th 
of March (all optical fluxes are host galaxy sub- 
tracted). This trend continued until June, where 
the source remained in a low optical state. Unfor- 
tunately with six nights the sampling of the opti- 
cal light curve during the MAGIC observations is 
quite low, hence an increase of the optical activity 
on a time scale of days can not be excluded. 

5.3. 1ES 1011+496 

The standard analysis performed for the whole 
source sample yields a significance of 2.5 a, which 
gives already a hint for a possible signal. In a more 
refined analysis, the cut in 2 , which determines 
the signal region, was reduced to 8 = 0.20°. In 
case of a weak signal the increased signal to back- 
ground ratio would lead to a higher significance. 
Also the SIZE-dependent cut for the background 
suppression is changed to a slightly higher value. 
The same analysis, performed on a data sample 
of the Crab Nebula, results in almost the same 
significance as with the standard coefficients, but 
with a ~ 13% lower 7-rate due to the reduced 7 
acceptance and a 37% lower background rate. 

This analysis yields 3.5 a for 1ES 1011+496 
which - if interpreted as a detection - corre- 
sponds to an integral flux of F(> 180 GeV) = 
(1.26 ± 0.40) x 10" 11 photonscm" 2 s" 1 . Further 
observations with the MAGIC telescope, triggered 
by an optical outburst in March 2007, shows a 
clear signal of 6.2 a within 18.7 hrs of observa- 
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Table 3 
Results of the analysis. 







Exp. 


Zd a 


^thrcs 










Sig. 


UL 


UL 


Source 


Mode 


[hrs] 


[degr.] 


[GeV] 


Excess c 


Backgr. 


Scale 


a 


c.u. d 


f.u. c 


1ES 0120+340 


w 


14.9 


12.2 


190 


-48 


5358 





.33 


-0.6 


0.032 


0.75 


RX J0319. 8+1845 


w 


4.7 


14.3 


190 


9 


2225 





.33 


0.2 


0.049 


1.15 


RX J0319. 8+1845 


on 


6.5 


14.2 


190 


-95 


3257 





.86 


-1.2 


0.033 


0.78 


1ES 0323+022 


w 


11.4 


29.0 


230 


55 


5262 





.33 


0.7 


0.064 


1.16 


1ES 0414+009 


w 


17.8 


29.7 


230 


176 


7309 





.33 


1.8 


0.057 


1.03 


1ES 0806+524 


w 


17.5 


26.8 


230 


111 


6174 





.33 


1.2 


0.056 


1.01 


1ES 0927+500 


w 


16.1 


22.1 


230 


72 


5721 





.33 


0.8 


0.052 


0.94 


1ES 1011+496 


w 


14.5 


23.6 


230 


200 


4857 





.33 


2.5 


0.086 


1.55 


1ES 1218+304 


w 


14.6 


26.6 


230 


400 


5423 





.33 


4.6 


0.073 f 


1.31 f 


RX J1417.9+2543 


on 


13.0 


9.7 


190 


-137 


9007 


1 


.03 


-1.0 


0.023 


0.54 


1ES 1426+428 


on 


6.1 


16.6 


190 


-7 


2561 





.24 


-0.1 


0.050 


1.18 


RX J1725. 0+1152 


on 


5.3 


17.4 


190 


-69 


2001 





.98 


-1.1 


0.046 


1.08 



a Mean zenith angle of the observation 

b Peak response energy for a Crab like spectrum 

c Background subtracted signal events for 6 < 0.23° 

d Integral flux above the threshold given in units of the flux of the Crab Nebula (crab units, c.u.) 
°Integral flux above the threshold given in flux units (f.u.) of 10~ n photons cm -2 s _1 
integral flux above threshold in cu. and f.u. corresponding to the 4.6 a excess 



8 



tion, resulting in flu x ~ 50% higher than in 2006 
([Albert et al.ll2007ef l. 



5.4. 1ES 1426+428 

The VERITAS collaboration reported a steep 
spectrum above 300 GeV for their observations in 
2001, well fitted b y a power law wit h spectral in- 
dex -(3. 50 ±0.35) (|PeteieEaD[2003). Extrapolat- 
ing the spectral fit to 200 GeV, it yields an integral 
flux of 0.50 Crab above 200 GeV, which is by a fac- 
tor of 10 larger than the u.l. presented in this work. 
Previous measurements yield a marginal detection 
in 2000 and u pper limits for the data taken from 
1995 to 1999 (jHoran et alj |2002| ) with the most 
stringent one of 0.08 Crab above 350 GeV. 

The HEGRA collaboration published a much 
harder spectrum at higher energies (above ~ 800 GeV) 
for their combined 1999 and 2000 data, which 
were well fitted b y a power law with sp ectral in- 
dex -(2.6 ± 0.6) (jAharonian et alj|2002h . An ex- 
trapolation of the power law yields an integral flux 
above 200 GeV of 0.075 Crab. Due to the large ex- 
trapolated energy range, combined with the large 
statistical error of 0.6 for the slope, the uncertainty 
is a factor of two. Further measurements in 2002 
with the HEGRA telescopes sh owed the source 



in a 2.5 times lower flux state ( Aharonian et al 



2003) 



The u.l. for the flux above 200 GeV presented 
in this work, indicates a lower flux than measured 
from 1999 to 2001 during several campaigns with 
different telescopes, whereas it is consistent with 
the low flux level observed in 2002. 

6. Spectral Energy Distribution 

The complete set of sources, described in 
Sect. 2, amounts to 14 objects (without 1ES 
1727+502) and includes the five established TeV 
sources Mrk 501, 1ES 1218+304, 1ES 1426+428, 
Mrk421 and 1ES 2344+514. To better under- 
stand the spectral energy distribution of VHE 
gamma ray emitting blazars, the sample is en- 
larged by including ten more sources which fullfill 
the same selection criteria in X-ray flux and red- 
shift, but which deviate only in the zenith angle 
cut. This means, they can only be observed under 
large zenith distances from the MAGIC site or 
not at all, in which case there is still information 
available from corresponding H.E.S.S. observa- 



tions. An exception is PG 1553+113 where the 
redshift is not known. For the total, enlarged 
sample of 24 HBLs, multi-wavelength data are 
collected in the foll owing bands: radio (5 GHz, 



Donato et al. ( 200ll )l, optical (R-band, 640nm, 
simultaneous data from KVA or, if not avail- 



able from iDonato et all (|200lh ). X-rays (IkeV, 
iDonato et all (|200ll )). and 7-rays (200 GeV, see 
TablelH TableGD and ref. in Table[T]). For 1ES 
1426+428, in addition to the upper limit derived 
in this work, the extrapolation of the spectrum, 
as measured by HEGRA in 1999/2000, is used to 
mark the detected flux. The optical data were cor- 
rected for galactic extinction, using the coefficients 
from the NASA Extragalactic Dat abase (NED) 



which are calculated according to ISchlegel et al 
<|l998h . 



For the 7-ray flux, sizeable attenua tion is ex- 

pected from current model s of the MRF (jHauser fc Dwek 



2001b iKneiske et all 12004 ) . Therefore all u.l. at 



200 GeV as well as the measured fluxes of the 
detected HBLs are corrected for the absorption 
by multiplying with exp(r(200 GeV, z)) where r 
denotes the pair production optical depth. The 
"best-fit 2006" MRF-model (Kneiske 2007, in 
preparation) is employed. This model is based o n 
the "best fit" model from IKneiske et all (|2004 ). 
but with a lower star formation rate to keep 
the energy density in the optical band closer to 
the lower limits, derived from the galaxy number 
co unts. It is also consistent with the u.l. derived 
by [Ah aronian et all (|2006al ) from the VHE spec- 
trum of 1ES 1101-233 (z = 0.186). The optical 
depth values for 200 GeV photons for all sources 
of this sample are listed in TableO 

All fluxes are K-corrected. Radio spect ral in- 
dicejl of ten of the sources can be found in Landt 
(|2003l ). For the other 14 objects the average 
value cur = 0.23 of the ten sources is used. 
For the optical data, the spectral indices of nine 
sources, calcula ted at slightly higher w avelengths, 



are taken from iBersanelli et all (|1992f ). For the 
other 15 objects, the average value ao = 0.65 of 
the nine sources is use d. At 1 keV, t he sp ectral 
indices are taken from IDonato et all (|200ll ). ex- 
cept for 1ES 0229+200, which is not included in 
this compilation. Instead, the flux is taken from 



2 Spectral index defined by the photon flux F oc v a [pho- 
tons cm~ 2 s" 1 ! 
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Costamante fc Ghisellini (2002) together with the 
average value for the spectral index ax = 1-36 
of all other sources. At 200 GeV the measured 
spectral indices are used for the detected sources, 
while for the non-detected ones the average value 
a 7 = 2.0 is used. To take into account the en- 
ergy dependent attenuation at VHE, which causes 
a hardening of the spectrum, the measured spec- 
tral indices are changed by -0.4 for 0.1 < z < 0.2, 
-0.8 for 0.2 < z < 0.3 and remain unchanged for 
z < 0.1. 

A special treatment is necessary to derive the 
flux at 200 GeV from 1E S 0229+200, recently dis - 



covered at VHE 7-rays (jAharonian et a l. 2007c) 



The spectrum is measured above 580 GeV, well fit- 
ted by a power law with spectral index —(2.51 ± 
0.19). As the source is located at z — 0.1396, 
strong absorption is expected at these energies. 
Therefore the spectrum is first deabsorbed and 
afterwards extrapolated to lower energies (Ta- 
blc[4]). The resulting intrinsic spectrum is well 
described by a power law with a spectral index 
of -(1.09 ± 0.25) (flux normalisation: (4.24 ± 
0.81) x 10~ 12 cm- 2 s- 1 TeV" 1 at ITeV). This re- 
sult is in good agreem ent with the results from 
IStecker fc Scully ( 2007 ). yielding model depen- 
dent intrinsic spectral indices in the range from 
1.1 ±0.3 to 1.5 ±0.3. 

After these corrections, the broad-band spec- 
tral indices^ &1-2 between the different energy 
regimes are calculated. Also the luminosities vh v 
are calculated assuming isotropic emission and 
with the use of the following cosmological param- 



eters: 
SZ™ — 



#0 

0.27. 



71kms- 1 Mpc~ 



0.73, and 



Discussion 



7.1. Gamma ray emitting HBLs? 

One may ask whether the gamma ray emitting 
HBLs can be distinguished from other HBLs based 
on their spectral energy distributions. Finding the 
answer is hampered by a number of problems. The 
variable peak frequencies, which cannot be easily 
detected in fixed energy bands, the gamma ray 
attenuation due to pair production in the meta- 
galact ic radiation field, a nd the flux variability. 
From iDonato et al.l (|200ll ). the amplitude of the 



-log(Fi/F 2 )/log(!/i/^ 2 ) , v\<V2 



flux variability at 1 keV amounts to a factor of six 
for the sources with multiple entries in the cat- 
alogue. Similar or larger amplitudes can be ex- 
pected at VHE. However, our sample is not trig- 
gered by flux variability, the duty cycle of flares 
generally seems to be rather low, and the observed 
fluxes or flux upper limits may therefore be char- 
acteristic of the quiescent average fluxes. 

FigurcU] shows the broad-band spectral index 
ctRO vs. ctox for all 24 HBLs as described in the 
previous section. The distribution is quite ho- 
mogeneous. As the data are not simultaneously 
taken, the uncertainties due to flux variations have 
to be taken into account. In the case of aox a flux 
variability of a factor six at 1 keV coresponds to 
a change in the spectral index of 0.29 (if the op- 
tical flux remains the same). This is still below 
the difference of 0.6 between the lowest and high- 
est values of aox for the detected VHE sources. 
Thus there are significant differences in the SED 
among the HBLs studied here. As the variability 
in the radio and optical band for HBLs is lower 
than at X-rays or VHE 7-rays and the change of 
the spectral index oro is smaller for different flux 
ratios, the variation of oro is much lower than for 
aox- However, it is not possible to distinguish be- 
tween sources detected at VHE and non-detected 
ones based on spectral indices. 

7.2. Gamma-to- X-ray spectral index 

To unveil the physical state of the emitting 
plasma, we seek to find the characteristic ratio 
of the two peaks in the spectral energy distri- 
bution of HBLs which is related to the ratio of 
photon and magnetic field energy densities. Fig- 
ure[5]shows the broad-band spectral index ao 7 vs. 
ax 7 - Both indices are distributed in a narrow 
band around unity for all detected sources, where 
a = 1 represents the case that the energy output 
in both frequency bands is the same. The con- 
straints from the u.l.s on the 7-ray flux can not ex- 
clude the region which is spanned by the detected 
TeV sources. In the framework of SSC models, 
the optical and X-ray band belong to synchrotron 
emission of highly relativistic electrons (first peak 
in the SED), while the VHE 7-rays are produced 
by inverse Compton scattering (second peak in the 
SED). If the magnetic energy density mb is equal 
to the photon energy density it p h, the energy out- 
put for the synchrotron and the Inverse Compton 
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Table 4 

hbls detected at vhe which do not belong to the sample described in sect. 2 . 



Source 


Z7 a 

try 


Kei. 


1ES 0229+200 


1.6 b 


Aharonian et al. (20070) 


1ES 0347-121 


4.25 


Aharonian et al. (2007b) 


PKS 0548-322 


1.9 


SuDerina et al. (2007) 


1ES 1101-232 


2.93 


Aharonian et al. (2006a) 


Mrk 180 


11.0 


Albert et al. (2006c) 


PG 1553+113 


11.5 


Albert et al. (2007a) 


1ES 1959+650 


17.4 


Albert et al. (2006a) 


PKS 2005-304 


6.63 


Aharonian et al. (2005a) 


PKS 2155-489 


26.3 


Aharonian et al. (2005b) 


H 2356-309 


2.78 


Aharonian et al. (2006a) 



a zAFV(200GeV) in units of 10~ 12 ergcm" 2 s" 1 

b This value comes from an extrapolation as de- 
scribed in the text 



emission as well as the peak luminosities are the 
same. In case of HBLs, the peak frequency of the 
synchrotron emission is always at higher frequen- 
cies than the optical band and below IkeV for 
most of them (except the extreme blazars with a 
hard spectrum at IkeV). At VHE energies, most 
of the detected sources show a soft spectrum above 
~ 200 GeV, indicating observed peak energies be- 
low 100 GeV. Due to absorption of VHE 7-rays 
in the MRF, the intrinsic peak energies could also 
reach several TeV for sources located at higher red- 
shifts. In that context, the scattering of ao 7 and 
ax 7 around unity could be explained by a con- 
tinuous distribution of peak frequencies for HBLs, 
measured with fixed bandwidth. 

Figurc[S] shows the broad-band spectral index 
ax-y vs. X-ray luminosity v-^Ly^. The average en- 
ergy output at 1 keV never exceeds significantly 
the one at 200 GeV {a Xl > 0.97). For half of the 
detected sources, the energy output in these bands 
is almost the same (ax-y — 1), while for the other 
ones the energy output at 200 GeV is significantly 
lower. There is a tendency that this effect shows 
up at higher X-ray luminosities. For four of the ten 
u.L, ax-y = 1 can not be excluded, whereas for the 
other six sources ax 7 has to be larger than unity. 
For the non-detected sources of the sample, fur- 
ther observations with longer exposures are needed 
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Fig. 4. — The broad-band spectral index ckro vs. 
ci^ox- The filled symbols mark the spectral in- 
dices of the sources which belong to the sample 
described in Sect. 2. They are further divided into 
detected (stars) and non-detected sources (trian- 
gles). The open circles represent all other HBLs 
detected at VHE 7-rays so far. 
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Table 5 

Upper limits on the 7-ray flux at 200 GeV under the assumption of a -3.0 power law 
spectrum togehter with the optical depth and the simultaneous optical data. 



J^thrcs 



Source 


Mode 


GeV 


F b 




z 




T C 




F d 






1ES 0120+340 


w 


170 


4.0 





272 





53 





47 + 





05 


10.5 


RX J0319.8+1845 


w 


170 


4.2 





190 





32 





48 + 





10 


4.3 


RX J0319.8+1845 


on 


170 


6.2 





190 





32 





14 + 





10 


4.3 


1ES 0323+022 


w 


190 


8.7 





147 





22 


1 


82 + 





19 


7.8 


1ES 0414+009 


w 


190 


7.7 





287 





57 










12.1 


1ES 0806+524 


w 


190 


7.6 





138 





21 


8 


00 + 





23 


11.9 


1ES 0927+500 


w 


170 


7.1 





188 





31 










9.7 


1ES 1011+496 


w 


170 


10.9/6. 5 f 





212 





37 


11 


49 + 





13 


5.2 


1ES 1218+3048 


on 


120 


10.1 





182 





29 


6 


13 + 





13 


21.2 


1ES 1218+304 


w 


190 


7.7 





182 





29 


4 


99 + 





11 


21.2 


RX J1417.9+2543 


on 


140 


2.5 





237 





43 


2 


11 + 





29 


8.7 


1ES 1426+428 


on 


140 


5.5 





129 





19 


1 


87 + 





15 


18.5 


RX J1725.0+1152 


on 


190 


6.3 


> 


0.17 





27 h 


13 


27 + 





09 


8.7 


1ES 2344+514 1 


w 


180 


11.5 





044 





05 


3 


37 + 





25 


12.0 



a Peak response energy for a power law spectrum with index -3.0 
b ^F„(200GeV) in units of 10~ 12 erg cm" 2 s" 1 
c optical depth t(z) at 200 GeV 

d i/-F , „(640 nm) in units of 10~ 12 erg cm -2 s _1 (Host galaxy substracted) 
c vF lJ {l'keY) in units of 10~ 12 erg cm~ 2 s _1 taken from lDonato et all ( 2001 ) 
f Flux under the assumption of a detection 



s Values from I Albert et all (|2006bl ) 

h For a redshift of 0.17 

'Values from I Albert et all (|2007bh 
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Fig. 5. — The broad-band spectral index ao 7 
vs. ax-y- The arrows mark the upper limits for 
the spectral indices, as calculated in this work, 
whereas the stars indicate the spectral indices of 
the detected sources that belong to the sample de- 
scribed in Sect. 2. The open circles mark all are 
other HBL objects detected a VHE 7-rays so far. 
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Fig. 6. — The broad-band spectral index ax-y vs. 
the X-ray luminosity vxLx- The arrows mark the 
upper limits for the spectral indices, as calculated 
in this work, while the stars indicate the values 
for the detected sources that belong to the sample 
described in Sect. 2. The open circles mark all 
other HBL objects detected at VHE 7-rays so far. 



to reach the axy = 1.12-line (corresponding to a 
nine times lower output at 200 GeV compared to 
1 keV) , which includes all HBLs detected at VHE 
so far. Note that if 0^7 = 1 f° r the peak fre- 
quencies is valid, the tendency of increasing ax-y 
values with increasing luminosities at 1 keV could 
be interpreted as a shift of the inverse Compton 
peak to lower values. Similar studies aiming at a 
ch aracteriza t ion o f VHE blazars are also pursued 
bv I Wagner" 



(120071) 



7.3. Constraints on the luminosity func- 
tion of 7-ray emitting HBLs 

Figure[7] shows the luminosity v 1 L 1 at 200 GeV 
vs. the redshift. All detected sources are above or 
within the line that marks the corresponding lumi- 
nosity to a flux of 2 jiiJy at 1 keV. The absorption 
of 7-rays by the MRF increases with redshift, so 
that at a redshift of z = 0.3, the emitted luminos- 
ity becomes twice larger than the measured one. 

The luminosity function at 7-ray energies of 
HBLs is poorly known, since there has not been 
a complete survey and the number of detected 
sources is still rather low. Nevertheless, we can 
try to constrain the VHE luminosity function from 



our observations. To this end we derive upper 
limits on the cumulative omnidirectional flux at 
200 GeV from X-ray bright HBLs below z = 0.3, 
noting that GLAST will measure the diffuse ex- 
tragalactic background up to 200 GeV. The se- 
lection criterion for the declination of the sam- 
ple corresponds the a patch of the sky with a 
size of 5.55 sr (or 44% of the sky). The sum 
over all 14 sources of the sample divided by the 
5.55 sr patch results in an upper limit on the to- 
tal intensity at 200 GeV of I V he(200 GeV) = 
£•2.76 x 10 _9 GeV cm -2 sr -1 s — 1 , where e accounts 
for the incompleteness of the sample. Assuming 
an isotropic distribution of HBLs, e should be 
larger than unity. We know about two sources, 
not included in this calculation (1ES 1727+502 
and 1ES 0229+20 0). There are also 5 HBLs 



from iDonato et al.l (|200lh that fullfill the criteria 
of declination and X-ray flux, but are located at 
higher redshifts. Including these sources the sam- 
ple would increase to 21 sources. With respect 
to the ~ 50% sky coverage of the Einstein Slew 
Survey, where most of the sources belong to, e = 3 
seems to be a reasonable assumption (factor 1.5 to 
take the "known" sources into account and a factor 
2 for the assumed sky coverage of 50% for the ex- 



13 



_r 



10* 

































o 


* 


V 






*^ 






























* 

c 


A 


Q 








: *> 
































0.05 



0.1 



0.15 

z 



0.2 



0.25 



0.3 



Fig. 7. — The 7-ray luminosity v 1 L~ j vs. redshift. 
The arrows mark the u.l.s calculated in this work, 
whereas the stars indicate the detected sources 
that belong to the sample described in Sect. 2. The 
triangle marks 1ES 1011+496 if interpreted as a 
detection. The open circles mark all other HBL 
objects detected at VHE 7-rays so far. The dashed 
line indicates a flux of 4.8x 10 -12 erg cm~ 2 s _1 cor- 
responding to 2 (itJy at 1 keV. The solid line corre- 
sponds to the same flux taking into account 7-ray 
attenuation at 200 GeV. 

tended sample of 21 sources). This would lead to 
an upper limit for the total intensity at 200 GeV of 
Ivhe(200 GeV) = 8.3 x lO^GeVcixr 2 sr" 1 s" 1 . 
Not that this u.l. is already conservative as also 
sources are included which are not in the Einstein 
Slew Survey (would result in a higher sky cover- 
age) and the u.l. is dominated by the two bright- 
est sources Mrk 501 and Mrk 421 (which would 
mean that three times more sources would not 
necess arily lead to a three times hig her flux). Re- 
cently iKneiske fc Mannheim ( 2007 ) showed that 
HBLs could contribute up to 30% of the extra- 
galactic background radiation at GeV energies, 
when including cascade emission from sources at 
higher redshifts. The luminosity function used for 
their calculation was derived from the X -ray lu- 
minosity function ( Beckmann et al. 2003f) assum- 



ing the same luminosity above 300 GeV as from 
0.5 keV to 2 keV. For the HBL contribution of faint 
point sources at 200 GeV, they obtain (somewhat 
model-dependent) ^km"^ (200 GeV) = (4 - 6) x 
10~ 8 GeV cm~ 2 sr -1 s _1 which is within the upper 



limit obtained here. For the total intensity, includ- 
ing the diffuse component due to electromagnetic 
cascading, their result is I^jf use '(200 GeV) = 
1.0 x 10~ 7 GeVcm~ 2 sr -1 s _1 . This estimate is 
based on assuming that the emitted VHE spec- 
tra generally have peaks at energies well in excess 
of 200 GeV. 

8. Conclusions 

Searching for VHE 7-ray emission from a sam- 
ple of twelve X-ray bright HBL objects, 1ES 
1218+30.4 at a redshift of z = 0.182 was dis- 



covered for the first time at VHE (Albert et al 



2006bl ). The already established VHE source 
1ES 2344 +514 has been det ected with high 
significance! Albert et al.ll2007rJ ). albeit in a state 
of low activity. For ten sources no significant sig- 
nal was seen, resulting in upper limits on their 
integral flux above ~ 200 GeV between 2.3% and 
8.6% of the Crab Nebula flux on a 99% confi- 
dence level. There is a hint for a signal from 
1ES 1011+496 on a 3.5 c level, which has been 
confirmed as a source of VHE 7-rays by a sec- 
ond MAGIC observ ation campaign, trig gered by 



a high optical state (jAlbert et al 



2007. 



With fixed-schedule observations a bias to flar- 
ing emission states was avoided, tacitly assuming 
that the duty cycle of the flares is short when com- 
pared with the exposure time. As shown in the 
case of 1ES 2344+514, quiescence does not nec- 
essarily preclude the sources from being detected. 
As a matter of fact, the upper limits obtained for 
the other sources still lie in a region of parameter 
space bracketed by the detected sources, and cor- 
respond to a VHE energy flux on the level of the 
X-ray energy flux with an indication of an increas- 
ing with increasing X-ray luminosity. It thus 
seems to be a question of time that more sensitive 
telescopes, and in particular those with a lowered 
energy threshold such as MAGIC-II, potentially 
capturing shifting peaks, will eventually lead to a 
detection of all known bright HBLs. It will be im- 
portant to obtain a complete catalogue of HBLs 
from the planned eROSITA and GLAST all-sky 
surveys to study them on a much larger statistical 
basis in the future. 

The detected sources deviate in no apparent 
pattern from the so- far non-detected sources. A 
spectral shape with two equal-height bumps is ex- 
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pected in Synchrotron Self Compton models for 
the case of balanced energy densities of photons 
and the magnetic field. The dispersion in the dis- 
tribution of the X-to-7-ray luminosity ratio (the 
lowest X-to-7-ray luminosity ratio of a detected 
HBL is 1/9) would then reflect variations of the 
peak position w.r.t. the observed band or the effect 
of flux variability (e.g., high state X-ray emission 
vs. quiescent VHE emission). 
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